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Abstract 

The applications of wireless power transfer technology to wireless communications can help build a wireless 
powered communication network (WPCN) with more reliable and sustainable power supply compared to the 
conventional battery-powered network. However, due to the fundamental differences in wireless information and 
power transmissions, many important aspects of conventional battery-powered wireless communication networks 
need to be redesigned for efficient operations of WPCNs. In this paper, we study the placement optimization of 
energy and information access points in WPCNs, where the wireless devices (WDs) harvest the radio frequency 
energy transferred by dedicated energy nodes (ENs) in the downlink, and use the harvested energy to transmit 
data to information access points (APs) in the uplink. In particular, we are interested in minimizing the network 
deployment cost with minimum number of ENs and APs by optimizing their locations, while satisfying the energy 
harvesting and communication performance requirements of the WDs. Specifically, we first study the minimum- 
cost placement problem when the ENs and APs are separately located, where an alternating optimization method 
is proposed to jointly optimize the locations of ENs and APs. Then, we study the placement optimization when 
each pair of EN and AP are co-located and integrated as a hybrid access point, and propose an efficient algorithm 
to solve this problem. Simulation results show that the proposed methods can effectively reduce the network 
deployment cost and yet guarantee the given performance requirements, which is a key consideration in the future 
applications of WPCNs. 
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I. Introduction 

Modem wireless eommunieation systems, e.g., eellular networks and wireless sensor networks (WSNs), 
are featured by larger bandwidth, higher data rate and lower eommunieation delays. The improvement on 
eommunieation quality and the inereased data proeessing eomplexity have imposed higher requirement on 
the quality of power supply to wireless deviees (WDs). Conventionally, WDs are powered by batteries, 
whieh have to be replaeed/reeharged manually onee the energy is depleted. Alternatively, the reeent 
advanee of radio frequeney (RF) enabled wireless power transfer (WPT) provides an attraetive solution 
to power WDs over the air [HI, [|2l. By leveraging the far-field radiative properties of mierowave, WDs 
ean harvest energy remotely from the RF signals radiated by the dedieated energy nodes (ENs) [|3|. 
Compared to the eonventional battery-powered methods, WPT ean save the eost due to manual battery 
replaeement/reeharging in many applieations, and also improve the network performanee by redueing 
energy outages of WDs. Currently, tens of mierowatts (/iW) RF power ean be effeetively transferred to 
a distanee of more than 10 meters O The energy is suffieient to power the aetivities of many low-power 
eommunieation deviees, sueh as sensors and RF identifieation (RFID) tags. In the future, we expeet more 
praetieal applieations of RF-enabled WPT to wireless eommunieations thanks to the rapid developments 
of many performanee enhaneing teehnologies, sueh as energy beamforming with multiple antennas ||4| 
and more effieient energy harvesting eireuit designs [|5|. 

In a wireless powered eommunieation network (WPCN), the operations of WDs, ineluding data trans¬ 
missions, are fully/partially powered by means of RF-enabled WPT ||6l- [[T4l . A TDMA (time division 
multiple aeeess) based protoeol for WPCN is first proposed in [|6l, where the WDs harvest RF energy 
broadeasted from a hybrid aeeess point (HAP) in the first time slot, and then use the harvested energy 
to transmit data baek to the HAP in the seeond time slot. Later, ||7| extends the single-antenna HAP in 
ii to a multi-antenna HAP that enables more effieient energy transmission via energy beamforming as 
well as more speetrally effieient SDMA (spaee division multiple aeeess) based information transmission 
as eompared to TDMA. To further improve the speetral effieieney, (HI eonsiders using full-duplex HAP 
in WPCNs, where a HAP ean transmit energy and reeeive user data simultaneously via advaneed self- 
interferenee eaneelation teehniques. Intuitively, using a HAP (or eo-loeated EN and information AP), 
instead of two separated EN and information aeeess point (AP), to provide information and energy aeeess 
is an eeonomie way to save deployment eost, and the energy and information transmissions in the network 

'Based on the product specifications on the website of Powercast Co. (http://www.powercastco.coml, with TX91501-3W power transmitter 
and P2110 Powerharvester receiver, the harvest RF power at a distance of 10 meters is about 40 /rW. 
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can also be more efficiently coordinated by the HAP. However, using HAP has an inherent drawback that 
it may lead to a severe “doubly-near-far” problem due to distance-dependent power loss That is, the 
far-away users quickly deplete their batteries because they harvest less energy in the downlink (DL) but 
consume more energy in the uplink (UL) for information transmission. To tackle this problem, separately 
located ENs and APs are considered to more flexibly balance the energy and information transmissions 
in WPCNs [|9ll- lfTT]l . In this paper, we consider the method using either co-located or separate EN and 
information AP to build a WPCN. 

Most of the existing studies on WPCNs focus on optimizing real-time resource allocation, e.g., transmit 
signal power, waveforms and time slot lengths, based on instantaneous channel state information (CSI, 
e.g., [|6l|-[l8l|). In this paper, we are interested in the long-term network performance optimization based 
mainly on the average channel gains. It is worth mentioning that network optimizations in the two 
different time scales are complementary to each other in practice. That is, we use long-term performance 
optimization methods for the initial stage of network planning and deployment, while using short-term 
optimization methods for real-time network operations after the deployment. Many current works on 
WPCNs use stochastic models to study the long-term performance because of the analytical tractability, 
especially when the WDs are mobile in location. Eor instance, [Q applies a stochastic geometry model 
in a cellular network to derive the expression of transmission outage probability of WDs as a function of 
the densities of ENs and information APs. Similar stochastic geometry technique is also applied to WPT- 
enabled cognitive radio network in [fTOll to optimize the transmit power and node density for maximum 
primary/secondary network throughput. However, in many application scenarios, the locations of the WDs 
are fixed, e.g., a sensor network with sensor (WD) locations predetermined by the sensed objects, or an 
loT (intemet-of-things) network with static WDs. In this case, a practical problem that directly relates 
to the long-term performance of WPCNs, e.g., sensor’s operating lifetime, is to determine the optimal 
locations of the ENs and APs. Nonetheless, this important node placement problem in WPCNs is still 
lacking of concrete studies. 

In conventional battery-powered wireless communication networks, node placement problem concerns 
the optimal locations of information APs only, which has been well investigated especially for wireless 
sensor networks using various geometric, graphical and optimization algorithms (see e.g., lfT5]l - lfT9]l ). 
However, there exist major differences between the node placement problems in battery-powered and 
WPT-enabled wireless communication networks. On one hand, a common design objective in battery- 
powered wireless networks is to minimize the highest transmit power consumption among the WDs to 
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satisfy their individual transmission requirements. However, sueh energy-eonservation oriented design is 
not neeessarily optimal for WPCNs, beeause high power eonsumption of any WD ean now be replenished 
by means of WPT via deploying an EN elose to the WD. On the other hand, unlike information 
transmission, WPT will not induce harmful co-channel interference to unintended receivers, but instead 
can boost their energy harvesting performance. These evident differences indicate that the node placement 
problem in battery-powered wireless communication networks should be revisited for WPCNs, to fully 
capture the advantages of WPT. 

In this paper, we study the node placement optimization problem in WPCNs, which aims to minimize 
the deployment cost on ENs and APs given that the energy harvesting and communication performances 
of all the WDs are satisfied. Our contributions are detailed below. 

1) We formulate the optimal node placement problem in WPCNs using either separated or co-located 
EN and AP. To simplify the analysis, we then transform the minimum-cost deployment problem 
into its equivalent form that optimizes the locations of fixed number of ENs and APs; 

2) The node placement optimization using separated EN and AP is highly non-convex and hard to 
solve. To tackle the non-convexity of the problem, we first propose an efficient cluster-based greedy 
algorithm to optimize the locations of ENs given fixed AP locations. Then, a trial-and-error based 
algorithm is proposed to optimize the locations of APs given fixed ENs locations. Based on the 
obtained results, we further propose an effective alternating method that jointly optimizes the EN 
and AP placements; 

3) Eor the node placement optimization using co-located EN and AP (or HAP), we extend the greedy 
EN placement method under fixed APs to solving the HAP placement optimization, which is 
achieved by incorporating additional considerations of dynamic WD-HAP associations during HAP 
placement. Specifically, a trial-and-error method is used to solve the WD-HAP association problem, 
which eventually leads to an efficient greedy HAP placement algorithm. 

Due to the non-convexity of the node placement problems in WPCNs, all the proposed algorithms are 
driven by the consideration of their applicabilities to large-size WPCNs, e.g., consisting of hundreds of 
WDs and EN/AP nodes. Specifically, we show that the proposed algorithms for either separated or co¬ 
located EN and AP placement are convergent and of low computational complexity. Besides, simulations 
validate our analysis and show that the proposed methods can effectively reduce the network deployment 
cost to guarantee the given performance requirements. The proposed algorithms may find their wide 
application in the future deployment of WPCNs, such as wireless sensor networks and loT networks. 
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Fig. 1. Schematics of a WPCN with separate ENs and APs. 


The rest of the paper is organized as follows. In Seetion II, we first introduce the system models 
of WPCN where the ENs and APs are either separated or co-located. Then, we formulate the optimal 
node placement problems for the two cases in Section III, and propose efficient algorithms to solve the 
problems in Sections IV and V, respectively. In Section VI, simulations are performed to evaluate the 
performance of the proposed node placement methods. Einally, we conclude the paper and discuss future 
working directions in Section VII. 


II. System Model 

A. Separated ENs and APs 

Eor the case of separated ENs and APs, we consider in Pig. [T] a WPCN in consisting of M 
ENs, N APs and K WDs, whose locations are denoted by 2 x 1 coordinate vectors {uj|f = 1, • • • , M}, 
{vj|j = 1, • • • , A^}, and {wfc|/c = l,---,iC}, respectively. We assume that the energy and information 
transmissions are performed on orthogonal frequency bands without interfering with each other. Specifi¬ 
cally, the ENs are connected to stable power source and broadcast RE energy in the DL for the WDs to 
harvest the energy and store in their rechargeable batteries. At the same time, the WDs use the battery 
power to transmit information to the APs in the UL. The circuit structure of a WD to perform the above 
operations is also shown in Pig. [T1 

In a transmission block of length T, the M ENs transmit simultaneously on the same bandwidth in 
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the DL, where eaeh EN i transmits 

Xi{t) = t e[0,T], i = !,■■■ ,M. (1) 

Here, Pq denotes the transmit power, Si{t) denotes the pseudo-random energy signal used by the i-\h EN, 
whieh is assumed to be of unit power (E* [|si(t)p] = 1) and independent among the ENs (Ef [si(f)sj(t)] = 
0 if i 7 ^ j). The reason to use random signal instead of a single sinusoid tone is to avoid peak in transmit 
power speetrum density, for satisfying the equivalent isotropieally radiated power (EIRE) requirement 
enforeed by speetrum regulating authorities |[I1. Notice that the energy beamforming technique proposed 
in dUl is not used in our setup, as it requires accurate CSI and DE symbol-level synchronization, which 
may be costly to implement in a highly distributed WPCN network considered in this work. 
Accordingly, the received energy signal by the fc-th WD is 

yk{t) = k = (2) 

where ^ denotes the equivalent baseband channel coefficient from the i-th EN to the fc-th WD, which is 
assumed to be constant within a transmission block but may vary independently across different blocks. 
Eet hi^k — denote the channel power gain, which follows a general distribution with its mean 

determined by the distance between the EN and WD, i.e., 

E[hi^k] =/3\\ui-Wk\\~'^°, i = 1,-■ ■ ,M, k = I,-■ ■ K, (3) 

where do > 2 denotes the path loss exponent in DE, || ■ || denotes the / 2 -norm operator, and (3 = 

( 8 \ 

j with Ad and fd denoting the downlink antenna power gain and carrier frequency, respectively 
[l20]| . Then, each WD k can harvest an average amount of energy from the energy transmission within 
each block given by [O 

Qk=vTE[\yk{t)\^\{hi,k}] =vTPo(zf=ihi,k), k = l,---,K, (4) 

where 77 G ( 0 , 1 ] denotes the energy harvesting circuit efficiency, and the expectation is taken over 
the pseudo-random energy signal variations under fixed /ij^fc’s within the transmission block. Eet = 
E [Qk] /T denote the average energy harvesting rate over the variation of wireless channels (hi^s) in 
different transmission blocks, we have 

Afc = ?7/3 Po ■ E»=illw - k = l, ■■■,!<:. (5) 

In the UE information transmissions, we assume that each WD transmits data to only one of the APs. 
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To make the placement problem tractable, the WD-AP associations are assumed to be fixed, where each 
WD k transmits to its nearest AP regardless of the instantaneous CSI, i.e., 

jk = arg min | |vj - w^l |, /c = 1 • ■ ■ , A'. (6) 

7=1,■■■A 

Here, we assume no co-channel interference for the received user signals from different WDs, e.g., the 
WDs transmit on orthogonal channels. Besides, for the simplicity of analysis, we assume no limit on the 
maximum number of WDs that an AP could receive data from. Then, the average power consumption 
rate for WD k is modeled as 

/ifc = ^ ai,fc + a2,fc||vj^ - = (7) 

where Oi ^ denotes the constant circuit power of WD k, denotes the average transmit power as a 
function of the distance between WD k and its associated AP jk, and du > 2 denotes the UL channel 
path loss exponent. Besides, a 2 ,k > 0 denotes a parameter related to the transmission strategy used in 
the UL communication!^ In general, the model in ([7]) indicates that the transmit power increases as a 
polynomial function of the distance between the transmitter and receiver to satisfy certain communication 
quality requirement, e.g., minimum data rate or maximum allowed outage probability, which is widely 
used for wireless network performance analysis IfTTI . If22]| . 

B. Co-located ENs and APs 

A special case of the WPCN that we consider in Fig. [U is when the ENs and APs are grouped into 
pairs and each pair of EN and AP are co-located and integrated as a hybrid access point (HAP), which 
corresponds to setting M = N and u* = Vj for i = 1, • • • , M. With the network model and HAP’s 
circuit structure shown in Pig. [21 a HAP transfers RE power in the DE and receives information in the 
UE simultaneously on different frequency bands. Although the use of HAPs is less flexible in placing 
the ENs and APs than with separated ENs and APs, the overall deployment cost is reduced, because the 
production and operation cost of a HAP is in general less than the sum-cost of two separate EN and AP. 

Por brevity, we reuse the notation u*, i = 1, • • • , M, to denote the location coordinates of the M HAPs. 
Given other parameters unchanged, the expression of the average energy harvesting rate of the k-\h 


^For example, given a receive signal power requirement F*, to achieve a target data rate and maximum allowed outage portability tpk for 
WD k, we have 02 ,fc = ( 3 ^ 0 ^ )when truncated channel inverse transmission II20I is used under Rayleigh fading 

channel, where Au and /„ denote the uplink antenna gain and carrier frequency, respectively, and Ei(a;) = Jj°° j ■ dt denotes the 
exponential integral function. 
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Fig. 2. Schematics of a WPCN with co-located ENs and APs (HAPs). 

WD is the same as that in ([5]). Meanwhile, the average power consumption rate fik can be obtained from 
(|7]) by replacing Vj with u* as follows. 

/ffc = ai,fc + a2,fc||u4 - fc = l, (8) 

where ik is the index of the HAP that WD k associates with, i.e., 

4 = arg min \\ui - Wk\\, k = 1, ■ ■ ■ , K. (9) 

,M 

C. Net Energy Harvesting Rate 

With the above definitions, the net energy harvesting rates of the WDs in both cases of separate and 
co-located EN and AP are given by 

u)k Afc /Xfc, k 1, • • • , K. (10) 

In practice, the net energy harvesting rate can directly translate to the performance of device operating 
lifetime (see e.g., [l23]|). Specifically, given an initial battery level C, the average time before the fc-th 
WD’s battery depletes is —Cjuk when Uk < 0, and +00 when Uk > OO In other words, given a minimum 

^We neglect in this paper the battery degradation effect 1241 caused by repeated charging and discharging. 


December 10, 2015 


DRAFT 

















device operating lifetime requirement > 0, it must satisfy Uk > —C/T^ if < oo, and ccfc > 0 if 
Tk = oo. 


III. Problem Formulation 

In this paper, we assume that the locations of the WDs are known and study the optimal placement 
of ENs and information APs, which are either separated or co-located in their locations. This may 
correspond to a sensor network with sensor (WD) locations predetermined by the sensed objects, or an 
loT network with static WDs. In particular, we are interested in minimizing the deployment cost given 
that the net energy harvesting rates of all the WDs are larger than a common prescribed value 7, i.e., 
> 7 , k = 1 ■ ■ ■ , K, where 7 is set to achieve a desired device operating lifetime. 

A. Separated ENs and APs 

When the ENs and APs are separated, the total deployment cost is ciM + C 2 N if M ENs and N APs 
are used, where ci and C 2 are the monetary costs of deploying an EN and an AP, respectively. To solve 


the minimum-cost deployment problem, let us first consider the following feasibility problem: 

Find = [ui, • • • ,um] , = [vi, • • ■ , vjv] (11a) 

s.t. Afc (U^^) -/ifc > 7, k = l,---,K, (11b) 

< Uj < b^, i = 1, - ■ ■ , M, (1 Ic) 

b'<v,<b^ j = (lid) 

where Xk and pk are functions of and given in ([5]) and (|7]), respectively. The inequalities in (II Icl) 
and (II Idl) denote element-wise relations. Besides, |b^, b^} specifies a feasible deployment area for both 
the ENs and APs in which is large enough to contain all the WDs, i.e., b* < < b^, k = 1 ■ ■ ■ , K. 

Evidently, if (fTTI) can be efficiently solved for any M and N, then the optimal node placement solution 


to the considered minimum-cost deployment problem can be easily obtained through a simple two- 
dimension search over the values of M and N, i.e., finding a pair of feasible (M, N) that produces the 
lowest deployment cost ciM + C 2 N. 

Eor a pair of fixed M and N (M > 0 and > 0), we can see that (fTTI) is feasible if and only if the 
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the following 

problem is 

no smaller than 7, i.e.. 




max 

JJM 


min {Xk (U^) - 


77} 

( 12 a) 

s. t. 


b* < Uj < b^, i = 

= 1 , 

••• ,M, 

( 12 b) 



.c " 

VI 

> 

VI 

= 1 , 

••• ,iV. 

( 12 c) 

as its equivalent epigraphic form |l27l|, i.e.. 




max 

t 











s. t. 

Xk (u^) 

- dk >t, fc = 1, • 

’ * * 5 

K, 

(13) 


b* < Uj < b^, i = 1, • • • , M, 

b' < Vj < b'*, j = !,■■■ ,N. 


Given fixed M and N, (fTTI) is feasible if and only if the optimal objeetive of (fT3l) satisfies t* > 7. Then, 
the key difficulty of solving the optimal deployment problem is to find efficient solution for problem 

dUl). 


B. Co-located ENs and APs 

When the ENs and APs are integrated as HAPs, the total deployment cost is C 3 M if M HAPs are 
used. Here, C 3 denotes the cost of deploying a HAP, where in general C 3 < ci + C 2 . Similar to the case 
of separated ENs and APs, the minimum-cost placement problem can be equivalently formulated as the 
following feasibility problem for any fixed number of M > 0 HAPs, 

bind = [ui, • • • ,um] 

s.t. - /ifc (U"") > 7 , k = lr--,K, (14) 

b^ < Uj < b^, i = 1 , • • • , M, 

where Xk (U^) and /i^ (U^) are given in ([5]) and ([ 8 ]), respectively. Notice that the study on co-located 
ENs and APs is not a special case of that of separated ENs and APs. In fact, it adds extra constraints 
(Uj = Vj, 1 = 1 , • • • , M) to (fTTI) . which leave less flexibility to the nodes placement design and make 
the problem more challenging to solve. Equivalently, the feasibility of (fT4l) can be determined by solving 
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the following optimization problem 

max t 

t,u*f 

s.t. k = l,---,K, (15) 

< Uj < b^, i = 1 , - ■ ■ , M, 

and then eomparing the optimal objeetive t* with 7, to see whether > 7 holds. In the following 
Seetions IV and V, we propose effieient algorithms to solve problems (fTSl) and (fTSl) . respeetively. It is 
worth mentioning that the placement solution to (fT3l) and (fT5l) can be at arbitrary locations. When an 
EN (or a HAP) is placed at a location very close to an WD, the far-field channel model in ([3]) may 
be inaccurate. However, we learn from (fTSl) and (fTSl) that the optimal value t* is determined by the 
performance-bottleneck WD that is far away from the ENs and APs (i.e., the channel model in (3) 
applies practically), thus having very low energy harvesting rate and high transmit power consumption. 
Therefore, the potential inaccuracy of ([3]) will not affect the objective values of (fTSI) and (fTSl) . and the 
proposed algorithms in this paper are valid in practice. 

IV. Placement Optimization of Separated ENs and APs 

In this section, we study the node placement optimization for separately located EN and AP in problem 
(fTSl) . Specifically, we first study in Section IV.A the method to optimize EN placement assuming that 
the locations of APs are fixed in a WPCN. In Section IV.B, we further study the method to optimize the 
placement of APs given fixed EN locations. Based on the obtained results, we then propose in Section 
I VC an alternating method to jointly optimize the placements of ENs and APs. In addition, an alternative 
local searching method is considered in Section IV.D for performance comparison. 

A. EN Placement Optimization with Fixed AP Location 

We first consider the optimal EN placement problem when the locations of the APs are fixed, i.e., 
v/s are known. In this case, the WD-AP association jk is known for each WD k from ®, and /i^’s can 
be calculated accordingly from dV]). It is worth mentioning that the proposed algorithms under the fixed 
AP setup can be directly extended to solve EN placement problem in other wireless powered networks 
not necessarily for communication purpose, e.g., a sensor network whose energy is mainly consumed on 
sensing and processing data, as long as the energy consumption rates /i^’s are known parameters. With 
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v/s and being fixed, we can rewrite (fT3l) as 

max t 

t,u*f 

s-t- ■ E*=il|ui - -/ifc > t, k = l,---,K, (16) 

< Uj < b^, i = 1, • • • , M, 

where = ril3Po. We can see that (fT^ is a non-convex optimization problem, because ||uj — is 

neither a convex nor concave function in u^. As it currently lacks of effective method to convert (fT^ into 
a convex optimization problem, the optimal solution is in general hard to obtain. However, for a special 
case with M = 1, i.e., placing only one EN, the optimal solution is obtained in the following. By setting 
M = 1, (fT^ can be rewritten as 

max t (IVa) 

t,Ui 

s. t. ||ui - = (17b) 

t + Hk 

< ui < b^. (17c) 

Although (fTTI) is still a non-convex optimization problem (as + Hk) is not a concave function in t), 
it is indeed a convex feasibility problem over ui when t is fixed, which can be efficiently solved using 
the interior point method (ITtII . Therefore, the optimal solution of (fTTI) can be obtained using a bi-section 
search method over t, whose pseudo-code is given in Algorithm [TJ Notice that the right hand side (RHS) 
of (llVbl) is always positive during the bisection search over t G {LBi,UBi). Besides, we can infer 
that Algorithm [T| converges to the optimal solution t*, because problem (fTTI) is feasible for t < t* and 
infeasible otherwise. The total number of feasibility tests performed is log 2 [{UBi — LBi) /cxi], where 
(Ti is a predetermined parameter corresponding to a solution precision requirement. 

Since placing one EN optimally is solved, we have the potential to decouple the difficult EN placement 
problem (fT^ into M relatively easy problems with M > 1. This motivates a greedy algorithm, which 
places the ENs iteratively one-by-one into the network. Intuitively, an optimal deployment solution of 
(fT^ should “spread” the M ENs among the K WDs to maximize the minimum energy harvesting rate. 
However, the optimal solution obtained from solving (fTTI) tends to place the single EN around the center 
of the cluster formed by the K WDs. Inspired by the optimal solution structures of (fT^ and (fTTI) . we 
propose in the following a cluster-based greedy EN placement method, where the newly placed EN 
optimizes the net energy harvesting rates of an expanding cluster of WDs, until all the WDs are included. 
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Algorithm 1: Bi-section search for single EN placement. 

input : WD locations w^’s, power consumption rates fik’s 
output: the optimal location of the EN u* 

1 Initialize: LBi < -maxfc^i,... /rfc, UB\ ■«— Po; 

2 repeat 

3 t ^ (t/Bi-hLPi)/2; 

4 if Problem dlZI l is feasible given t then 

5 LB\ •<— t\ 

6 ui ^ a feasible solution of GI1> given t\ 

7 else 

8 I [/Pi 3— f; 

9 end 

10 until \UBi — LB\\ < ai; 

11 Return u* <— the last feasible solution ui; 


In practice, we geographically partition the K WDs into M non-overlapping clusters (assuming K > 
M), denoted by {Wj|/ = 1, • • • , M}o This can be efficiently achieved by, e.g., the well-known iT-means 
clustering algorithm [|25]| . Then, in the Tth iteration (/ > 1), we obtain the optimal location of the /-th 
EN, denoted by u*, by maximizing the net energy harvesting rates of the WDs in the first i clusters as 
follows. 


max 

ii,Ui 

ti 

(18a) 

S. t. 

{U + fik- Xi-i,k) • llw - < ip, 



e {Wi u • • ■ u Wi}, 

(18b) 


b' < Ui < b\ 

(18c) 


where Xi-i,k denotes the accumulative RE power harvested at the fc-th WD due to the {i — 1) previously 
deployed ENs, given by 


A 


i—l^k 


0 , 


i = 1, 


(19) 


• Ep\||u* - Wfcll / > 1. 

In each iteration, u* can be efficiently obtained using a bi-section search method over similar to that 
in solving (fTTl) . The pseudo-code of the greedy algorithm is given in Algorithm [2l Notice in line 8 of 
Algorithm [2l the corresponding inequality in (I18bl) holds automatically when tj -f /ifc — Qi-i,k < 0 for 
some k, given a fixed G {LB 2 ^UB 2 ). Therefore, the corresponding constraints can be safely ignored 
without affecting the feasibility of (fTSl) . 


^The proposed node placement algorithms can apply to any clustering method used. Besides, the algorithm complexity is not related to 
the clustering method as long as the partitions of the WDs are given. For simplicity, we use if-means clustering algorithm in this paper. 
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5 


WDs not considered in the current EN deployment 
WDs considered in the current EN deployment 


g 


the EN currently being deployed 
ENs already deployed 


Fig. 3. Illustration of greedy algorithm for placing M = 3 ENs (with fixed APs). 


The greedy EN plaeement method is illustrated in Fig. [3l In this example, we first divide the WDs 
into M = 3 elusters in Fig. [3ta), then plaee the 3 ENs one-by-one in Fig. [3tb)-(d). When plaeing the 
1st EN or ENi, the algorithm only eonsiders the reeeived energy of the WDs in the 1st eluster (shaded 
WDs in Fig. Sb)) from the EN to be placed; for the 2nd EN or EN 2 , it considers the received energy 
of the WDs in the 1st and 2nd clusters from the first 2 ENs; for the last EN or EN 3 , it considers the 
received energy of all the WDs from the 3 ENs. Notice that our greedy algorithm allows multiple ENs 
to be placed in the same cluster, because the placement of the i-th EN considers all the WDs in the first 
i clusters. 

Algorithm [21 applies Algorithm [T] M times, one for placing each EN, thus the total number of 
feasibility tests performed is Mlog 2 [{UB 2 — LB 2 ) /cr 2 ], where (T 2 is a parameter corresponding to a 
solution precision requirement. Besides, the time complexity of solving each convex feasibility test using 
the interior point method is O [\/K + 2Mlog {K + 2M)) [|27]| . Therefore, the overall time complexity 
of Algorithm |2] is (9 [My/K + 2M log {K + 2M)), which is moderate even for a large-size network 
consisting of, e.g., tens of ENs and hundreds of WDs. 
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Algorithm 2: Greedy algorithm for M EN placement. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


input : WD locations w^’s, N AP locations v/s; 
output: locations of M ENs {uj, • • • , u^}; 

Initialization: Clustering the WDs into {Wi, i — 1, - • ■ , M} ; 

With v/s, calculate using and 0; 

for i = 1 to M do 

LB2 i -■ maxfc^i_... UB2 <— MPo ; 

Update {Qi-i^k, fc = 1 ■ ■ ■ ,K} using 

repeat 

ti ^ {UB2 + LB2)/2- 

Ignore the constraints in dl8bt with ti + Hk — Qi-i,k < 0; 
if Problem l UiSI I is feasible given U then 
LB 2 ^— ti\ 

u* a feasible solution of dlSb : 

else 

I UB 2 <— ti\ 

end 

until \UB2 — I/B 2 I < 0 - 2 ; 

end 

Return {u*, • • • , u^} 


B. AP Placement Optimization under Fixed ENs 

We then study in this subsection the method to optimize the placement of APs given fixed EN locations, 
i.e., Uj’s are known. In this case, A^’s are fixed and can be calculated using ([5]). With A^’s being fixed 
parameters, we can substitute © into (fT3l) and formulate the optimal AP placement problem under fixed 
ENs as follows 

max t 

s. t. Xk-ai^k - a 2 ,k\\vj^ - Wk\f^ >t, k = 1,-■ ■ ( 20 ) 

where jk is the index of AP that WD k associates with given in ®. The above problem is non-convex 
because of the combinatorial nature of WD-AP associations, i.e., jk's are discrete indicators. However, 
notice that if j^’s are known, (l20l) is a convex problem that is easily solvable. In practice, however, j^’s 
are revealed only after (l20l) is solved and the placement of APs is obtained. To resolve this conflict, we 
propose in the following a trial-and-error method to find feasible j^’s and accordingly a feasible AP 
placement solution to (l20l) . The pseudo-code of the method to solve (l20l) is presented in Algorithm [3] and 
explained as follows. 

As its name suggests, we first convert (l20l) into a convex problem by assuming a set of WD-AP 
associations, denoted by k = 1, - ■ ■ ,K, and then solve (l20l) for the optimal AP placement based on 
the assumed Next, we compare with the actual WD-AP associations after the optimal AP 
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placement is obtained using ®, denoted by k = 1, ■ ■ ■ , K. Speeifieally, we eheck if 
If yes, we have obtained a feasible solution to (l20l) : otherwise, we update k = 1, ■ ■ ■ , K 

and repeat the above process until = j^\ V/c. The eonvergenee of Algorithm [3] is proved in the 
Appendix and the eonvergenee rate is evaluated numerieally in Fig. of Seetion VI. Intuitively, the 
trial-and-error method is convergent because the optimal value of (l20l) is bounded, while by updating 
= jl°'\ we ean always improve the optimal objective value of (l20l) in the next round of solving it. 
As we will show later in Fig. of Seetion VI, the number of iterations used until eonvergenee is of 
eonstant order, i.e., 0(1), regardless of the value of N or K. There, the time eomplexity of Algorithm 
[3]is O + 2A^log {K + 2N)), as it takes this time eomplexity for solving (l20l) in each iteration. 


Algorithm 3: Trial-and-error method for N AP placement 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


input : K WD locations w^’s, M EN locations Uj’s; 
output: locations of N APs, i.e., {vj, • • • , v^}; 

Initialization: 

Separate the WDs into N clusters, and place each AP at a cluster center. Use to denote the initial AP locations; 

With calculate using 

With Uj’s, calculate Afc’s using (0. Let flag -t— 1; 

while flag = 1 do 

Given ’s, solve ( 120b for optimal AP placement v* ’s; 

Given vj’s, calculate using ([^; 

if /®'‘ some k then 

I Update fc = I,--- , A"; 

else 

A local optimum is found, return {v* 
flag t- 0; 

end 

end 


C. Joint EN and AP Placement Optimization 

In this subsection, we further study the problem of joint EN and AP placement optimization. In this 
ease, we eonsider both the loeations of ENs and APs as variables, such that the joint EN-AP plaeement 
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problem in (fT3]) can be expressed as 


max t 


\r^M II II —dn 

— <d2,k 


,K, 

< Uj < b^, z = 1, • • 

■,M 

b' < V,- < b\ J = 1, • • 

■■,N. 


( 21 ) 


Evidently, the optimization problem is highly non-convex because of the non-convex function ||uj — 
Wfc| 1“'^° and the discrete variables jkS. Based on the results in Section IV.A and IV.B, we propose an alter¬ 
nating method in Algorithm |4] to solve (|2T]) for joint EN and AP placement solution. Specifically, starting 
with a feasible AP placement, we alternately apply Algorithms [2] and [3] to iteratively update the locations 
of ENs and APs, respectively. A point to notice is that Algorithm [2] (and Algorithm [3]) only produces a 
sub-optimal solution to (fT^ (and (l20l)l. thus the objective value of (l2T]) may decrease during the alternating 
iterations. To cope with this problem, we record the deployment solutions obtained in L > 1 iterations and 
select the one with the best performance. The impact of the parameter L to the algorithm performance is 
evaluated in Pig. [5] of Section VI. Given the complexities of Algorithms [2] and [3l we can easily infer that 
the time complexity of Algorithm |4] is O (LM\/iT -f 2Mlog {K + 2M) + L^/K + 2N log {K + 2N)). 


Algorithm 4: An alternating method for joint AP-EN placement. 

input : K WD locations w^’s, L iterations; 

output; Locations of M ENs and N APs, i.e., uLs and v^’s.; 

1 Initialize: Separate the WDs into N clusters, and place each AP at a cluster center. Use Vj ’s to denote the initial AP locations; 

2 for / = 1 to L do 

3 if Z is odd then 

4 I Given v^ ’s, solve HD for Ui’s using Algorithm 2; 

5 else 

6 I Given Ui’s, solve HD for Vj’s using Algorithm 3; 

7 end 

8 2i •<— minfc=i_... (Afc — ^.k), where \k and fik are in H} and H)> respectively.; 

9 u® •(— Ui, i = 1, - ■ ■ ,M, and vj*' ■«— vj, j = 1, - ■ ■ ,N. 

10 end 

11 p ^ argmax!=i,... zj; 

12 Return: u)^^’s and v^A’s. 

'' J 


D. Alternative Method 

Besides the proposed alternating method for solving (|2TI) . we also consider an alternative local searching 
method used as benchmark algorithm for performance comparison. The local searching algorithm starts 
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with a random deployment of the M ENs and N APs, i.e., u/s and v^’s, and cheeks if the minimum net 
energy harvesting rate among the WDs, i.e., 


7~) A • / sr^M I I 

Pr = mm ip ■ - Wfc 


I —do 


- ai,fc - 02,fc _mm 11 
3 1 r" 


\du 


k=l,-,K 

can be increased by making a random movement to u/s and v^’s that satisfy 


Ui,Vj,i = !,■■■ ,M, j = !,■■■ ,N, 


11 — 

A*=il|ui 


u, 


+ EL"- 


Vi-v^l 


< 03 


( 22 ) 


(23) 


where is a fixed positive parameter. If yes, it makes the move and repeats the random movement 
process. Otherwise, if Pr cannot be increased, the algorithm has reached a local maximum and returns the 
current placement solution. Several off-the-shelf local searching algorithms are available, where simulated 
annealing [l26ll is used in this paper. In particular, simulated annealing can improve the searching result 
by allowing the nodes to be moved to locations with decreased value of Pr to reduce the chance of being 
trapped at local maximums. Besides, we can improve the quality of deployment solution using different 
initial node placements, which are obtained either randomly or empirically, and select the resulted solution 
with the best performance. 


V. Placement Optimization oe Co-located ENs and APs 

In this section, we proceed to study the node placement optimization problem (fT5l) for the case of 
co-located ENs and APs. The problem is still non-convex due to which the optimal solution is hard to be 
obtained. Inspired by both Algorithms [2] and [3l we propose in this section an efficient greedy algorithm 
for HAP placement optimization. 


A. Greedy Algorithm Design 

The node placement optimization problem (fTSl) is highly non-convex, because the expression of problem 
(fTSl) involves non-convex function | |uj—w^l in (U^) and minimum operator over convex functions 

in pk (U^). Since its optimal solution is hard to obtain, a promising alternative is the greedy algorithm, 
which iteratively places a single HAP to the network at one time, similar to Algorithm [2] for solving 
(fT^ which optimizes the EN locations given fixed APs. However, by comparing problems (fTSl) and (fT^ . 
we can see that the algorithm design for solving (fT5l) is more complicated, because each pk is now a 
function of u/s, instead of constant parameter in (fT^ . 

Similar to the greedy algorithm in Section IV.A, we first separate the K WDs into M non-overlapping 
clusters, denoted by {yVi\i = !,■■■ , M}, and add to the network a HAP in each iteration. Specifically, 
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in the i-th iteration, given that the previous (i — 1) HAPs are fixed, we obtain the optimal loeation of 
the Ath HAP, denoted by u*, by maximizing the net energy harvesting rates of the WDs in the first i 
elusters. To simplify the notations, we also use Ai_i ^ as in Seetion IV.A to denote the aeeumulative RF 
harvesting power of the WD k from the previously placed (i — 1) HAPs, which can be calculated using 
(fT^ . Besides, let denote the energy consumption rate of the /c-th WD after the first (i — 1) HAPs 
have been placed, where 

{ +CX), i = 1, 

(24) 

ai,fc + a 2 ,fcminj=i,...,j_i ||u* - i > 1. 

Notice that the only difference between placing the Ath HAP and the Ath EN in Section IV.A is that /i^ 
is now a function of u* instead of a given constant. By substituting (l24l) into (fT^ . the optimal location 
of the Ath HAP is obtained by solving the following problem 


where 


max 

U 

(25a) 

s. t. 

{ti + • 1 Iw - Wfcl< ip, 



k e {WiU---UWi}, 

(25b) 


< Uj < u^, 

(25c) 


Pi^k = min (/ii_i,fc,ai,fc + a 2 ,k\\ui - Wfc||'^^) . 

(26) 


From (|2^ . we can see that a WD may change its association to the Ath HAP, if the newly placed HAP is 
closer to the WD than all the other {i — 1) HAPs that have been previously deployed. This combinatorial 
nature of WD-AP associations makes problem (1251) non-convex even if ti is fixed. In the following, we 
apply the similar trial-and-error technique as that in Section IV.B to obtain a feasible solution to problem 

(125]). 


B. Solution to Problem ( |231) 

The basic idea to obtain a feasible solution of (1251) is to convert it into a convex problem given ti, 
and then use simple bi-section search over U. The convexification of (1251) is achieved by a trial-and-error 
method similar to that used for finding feasible WD-AP associations proposed in Algorithm [S] That is, we 
iteratively make assumptions on WD-AP associations and update the optimal placement of the Ath HAP 
obtained from solving (1251) based on the assumptions in the current iteration. With a bit abuse of notations. 
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here we reuse u* in eaeh iteration as the optimal loeation of the i-th HAP given the eurrent WD-AP 
assoeiation assumptions. Speeifieally, we assume whether the WDs ehange their assoeiations after the Ath 
HAP is added, i.e., assuming either < Ui.fc + a2,fc||u- - or > ai^k + 02,^1 |u- 

for eaeh k. Then, given a fixed U, eaeh eonstraint on k in (I25bl) belongs to one of the following four 
eases: 

1) Case 1: If we assume that WD k does not ehange its WD-HAP association after the i-th HAP is 
placed into the WPCN, or equivalently < ai,fc + a2,A;||u* —we can replace the corresponding 
constraint in (I25bl) with 

(ti + - Aj_i,fc) • I |uj - Wfcl < (f. (27) 

With a fixed U, (l27l) is a convex constraint if C + — K-i,k > 0. 

2) Case 2: If we still assume < ai,k + a2,fc||u* — Wfc||'^^, while U + — Xi-i,k < 0 holds, 

we can safely drop the constraint in (I25bl) without changing the feasible region of Uj. 

3) Case 3.' On the other occasion, if we assume that WD k changes its WD-HAP association, or 

k'i-i,k > 0 , 1 ,k + 02 ,k\\^i — the corresponding constraint in (I25bl) becomes 

{U + ai,fc + a2,k\\ui - - Xi-i,k) • ||w - < ip, (28) 

which can be further expressed as 

llw - - — < 0. (29) 

02,k 02,k 

Notice that, given a fixed C, (l29l) is a convex constraint if + Oi ^ — Aj_i ^ > 0. 

4) Case A: Otherwise, if we assume > ai^k + a2,fc||u* ^ and U + ai^k — Xt-i^k < 0 

holds, (l29l) is a non-convex constraint, as the left-hand-side (LHS) of (l29l) is the difference of two convex 
functions. Nonetheless, we show that (|29|) can still be converted into a convex constraint in this case. Let 
us first consider a function 

z{x) = ~ (30) 

02,k 02,k 

where x > 0 and ti + Oi^k — Ai_i,fc < 0. We calculate the first order derivative of z{x) and find that z 
increases monotonically when 

(ti T Oi k Xi—i k') dn) 

^ ^ 1 _ : _ : _ _ _ 

02,k (du + do) 

and decreases monotonically if x < Ti^k- Notice that Ti^k > 0 and z(0) = —^ < 0 always hold. Therefore, 
although z(x) is not a convex function, z(x) < 0 can still be equivalently expressed as x < 6i^k, with 9i^k 


L/au 


'^i.k 1 


(31) 
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being some positive number satisfying z{6i^k) = 0. The value of 6^^ ^ can be efficiently obtained using 
many off-the-shelf numerical methods, such as the classic Newton’s method or bi-section search method. 
A close comparison between the LHS of (l29l) and z{x) in (1^ shows that, by letting x ^ ||uj — Wfc||, 
we can equivalently express (|2^ as a convex constraint 

||ui - Wfcll < (32) 

when ti + ai^k — K-i,k < 0 holds. 

To sum up, given a fixed U, we tackle the k-th constraint in (I25bl) using one of the following methods: 

1) Replace by ([271) if assuming < ai^k + a 2 ,k\\u* - and U + fii-i^k - ><i-i,k > 0; 

2) Drop the constraint if assuming < ai,k + a2,A;||u* — and U + — Xi-i^k < 0; 

3) Replace by (|29|) if assuming ^i-i^k > ct-i,k + ct-2,k\W*i ~ '^k\\^^ and -f ai^k — K-i,k > 0; 

4) Replace by (l32l) if assuming > ai^k + a 2 ,k\\u* — and U + ai^k — K-i,k < 0. 

After processing all the K constraints in (I25bl) . we can convert (|251) into a convex feasibility problem given 
a set of WD-HAP association assumptions and a fixed fj. Accordingly, the optimal placement of the i-th 
HAP (u*) under the assumptions, can be efficiently obtained from solving (1251) using a bi-section search 
method over U. Similar to the trial-and-error technique used in Algorithm [3l we check if the obtained 
u* satisfies all the assumptions made. If yes, we have obtained a feasible solution of (l25l) . Otherwise, 
we switch the violating assumptions, then follow the above constraint processing method to resolve (1251) 
for a new u*, and repeat the iterations until all the assumptions are satisfied. The above trial-and-error 
method converges. The proof follows the similar argument as given in the Appendix, which proves the 
convergence of the trial-and-error method used for solving problem (l20l) . Thus, this proof is omitted here. 

C. Overall Algorithm 

Since the placement of a single HAP can be obtained via solving (1251) . we can iteratively place the M 
HAPs into the WPCN. The pseudo-code of the revised greedy algorithm is presented in Algorithm [51 For 
example. Fig. |4] illustrates the detailed steps taken to place the 2nd HAP, or HAP2 (total 3 HAPs while 
the first HAP, or HAPi is already placed). Specifically, we first assume in Fig. Ilta) that all the WDs in 
the 1st cluster associate with HAPi, and the WDs in the 2nd cluster associate with HAP2, after HAP2 is 
added into the network. Then, we obtain in Fig. |5b) the optimal placement of the HAP2 based on the 
association assumptions made. However, the obtained location of HAP2 results in a contradiction with 
the association assumption made on WDi (assumed to be associated with HAPi). Therefore, we change 
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- Assumed WD-to-HAP association A HAP location before optimization 

('V') 

- Actual WD-to-HAP association A Optimal HAP location under association assumption 


Fig. 4. Illustration of greedy algorithm for placement optimization with co-located EN and AP (HAP). 


the association assumption of WDi to HAP 2 , and recalculate the optimal placement solution for HAP 2 
(Fig. Sc)). In Fig. I2d), the newly obtained location of HAP 2 satisfies all the association assumptions, 
thus the placement of HAP 2 is feasible. Following the similar argument in the Appendix, the association 
assumption update procedure converges, because the optimal objective value of (l25l) is non-decreasing 
upon each association assumptions update (lines 10 — 27 of Algorithmic]). After obtaining the location of 
HAP 2 , a feasible location of the 3rd HAP, HAP 3 , can also be obtained using the similar procedures as 
above. Besides, we can infer that the time complexity of Algorithm [5] is O {M\JK + 2Mlog {K + 2M)), 
because it places the M HAPs iteratively, while the trial-and-error method used to place each HAP needs 
O {^JK + 2Mlog {K + 2M)) complexity. 

VI. Simulation Results 

In this section, we use simulations to evaluate the performance of the proposed node placement methods. 
All the computations are executed by MATLAB on a computer with an Intel Core i5 2.90-GHz CPU and 
4 GB of memory. The carrier frequency is 915 MHz for both DL and UL transmissions operating on 
different bandwidths. In the DL energy transmission, we consider using Powercast TX91501-IW power 
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Algorithm 5: Greedy algorithm for HAP plaeement optimization. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


input : K WD locations w^’s 

output: locations of M HAPs {ui, • • • , u^} 

Cluster the WDs into {Wi, i = 1, - ■ • , M}; 
for 2 = 1 to M do 

for each WD k do 

Update Xi-i,k and using GD and (mi; 

Assume WD k satisfies condition (a) or (6): 

(а) < ai,k +a 2 ,i:||u* - 

(б) > ai.fc + a2.fc||u* - 

end 

StopFlag 0; 

repeat 

LB < - S, UB S, (5 is sufficiently large; 

repeat 

ti {UB + LB)/2; 

Given ti, convert l l25t into a convex problem using the procedures in Section V.B; 
if Problem 4251 1 is feasible given ti then 
I LB <— ti', u* ■<— a feasible solution of l l25t ; 
else 

I UB*^ti', 

end 

until \UB — LB\ < a, a is sufficiently smalt, 
if all the K assumptions are valid then 
I StopFlag •<— 1; the i-th HAP location •<— u*; 
else 

StopFlag <r- 0; 

For each WD violating the assumption, switch the assumption from (a) to (b), or (b) to (a); 

end 

until StopFlag — 1; 

end 

Return the HAP locations {u*, • ■ • , u^}. 


transmitter with Pq = IW (Watt) transmit power, and P2110 reeeiver with = 3 dB antenna gain 
and 7] = 0.51 energy harvesting effieieney. Besides, we assume the path loss exponent do = 2.2, thus 
(3 = 6.57 X 10“^. In the UL information transmission, we assume that du = 2.5, Oi ^ = 50/iW and 
02 ,fc = 1.4 ■ 10“® for fc = 1, • • • , K, where a 2 ,k is obtained assuming Rayleigh fading and the use of 
truneated ehannel inversion transmission iflOl with reeeiver signal power — 70dBm (equivalent to 18dB 
SNR target with -88dBm noise power) and outage probability of 5%. All the WDs, ENs and APs are 
plaeed within a 24m x 24m box region speeified by = (0, 0)^ and = (24, 24)^. Unless otherwise 
stated, eaeh point in the following figures is an average performanee of 20 random WD plaeements, eaeh 
with 77 = 60 WDs uniformly plaeed within the box region. 

A. Separated EN and AP Deployment 

We first evaluate the performanee of the proposed alternating optimization method (Algorithm H]) for 
plaeing separated ENs and APs. Without loss of generality, we eonsider N = 6 APs and show in Pig. 
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Fig. 5. Performance comparison of the separated AP and EN placement methods (N = 6). 


[5] the minimum net energy harvesting rate Pr in ([T?]) aehieved by Algorithm |4] when the loeations of 
APs are jointly optimized with those of different number of ENs (M). Evidently, a larger Pr indieates 
better system performanee. Eor the proposed alternating optimization algorithm (AltOpt), we show both 
the performanee with L = 10 and 20. Besides, we also eonsider the following benehmark plaeement 
methods 

• Cluster eenter method (CC): separate the WDs into M elusters and plaee an EN at eaeh of the eluster 
eenters. Similarly, separate the WDs into N elusters and plaee the N APs at the eluster eenters; 

• Optimize only EN loeations: the APs are plaeed at the N eluster eenters; while the EN placement 
is optimized based on the AP locations using Algorithm [21 

• Eocal searching algorithm (ES) method introduced in Section IV, where the initial EN and AP 
locations are set according to the CC method and the best-performing deployment solution obtained 
during the searching iterations is used. 

Evidently, we can see that the proposed alternating optimization has the best performance among the 
methods considered. Specifically, significant performance gain is observed for AltOpt over optimizing 
EN placement only. The ES method has relatively good performance compared to AltOpt, especially 
when M is small, but the performance gap increases with M due to the increasing probability of being 
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(a) Average no. of iterations tili convergence (K=60 WDs) 


4 



I 1 


O'-'-'-'-^^^^^-'- 

4 6 8 10 12 14 16 18 20 22 24 

No. Of APS (N) 

(b) Average no. of iterations tili convergence (N=8 APs) 


4 



O'-'-^^^- 

40 50 60 70 80 90 

No. of WDs(K) 


Fig. 6. The average number of WD-AP association assumptions made before convergence of Algorithm [3] as a function of (a) N under 
bxed K — 60; and (b) as a function of K under fixed N = 8. 


trapped at local maximums with a larger M. The CC scheme has the worst performance as it neglects the 
disparity of energy harvesting/consumption rates among the WDs and precludes the case where multiple 
AP/ENs can be placed in a cluster. In practice, Fig. [5] can be used to evaluate the deployment cost of each 
algorithm. For instance, when Pr > —O.lmW is required, we see that the AltOpt (L = 10) on average 
needs 6 FNs, the FS method needs 8 FNs, optimizing FN placement only requires 9 FNs, and the CC 
method needs 20 FNs, with the same number of information APs deployed (i.e., N = 6). The above 
results show that, when the FNs and APs are separated, significant performance gain can be obtained by 
jointly optimizing the placements of FNs and APs, especially for large-size WPCNs that need a large 
number of FNs and APs to be deployed. In addition, as optimizing only FN locations corresponds to 
a special case of Algorithm |4] with L = 1, we can see that the performance gain is significant when L 
increases from 1 to 10. However, the performance improvement becomes marginal as we further increase 
L from 10 to 20. In practice, good system performance can be obtained with relatively small number of 
alternating optimizations, e.g., L = 10 in our case. 

We then show in Fig. the convergence rate of Algorithm [3l for which the convergence is proved in 
an asymptotic sense in the Appendix. In particular, we plot the average number of iterations (WD-AP 
association assumptions) used until the algorithm converges. Here, we investigate the convergence rate 
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Fig. 7. Performance comparison of the co-located AP and EN (HAP) placement methods. 

when either the number of APs (N) or WDs (K) varies. With fixed K in Fig. |6ta), we see that the number 
of iterations used till eonvergenee does not vary signifieantly as N inereases. Similarly in Fig. [6tb), with 
a fixed N = 8, we do not observe significant increase of iterations when K increases from 40 to 90. 
Besides, all the simulations performed in Fig. use at most 7 iterations to converge. Therefore, we can 
safely estimate that the number of iterations used till convergence is of constant order, i.e., 0(1), which 
leads to the complexity analysis of Algorithm |4] in Section IV.C and Algorithm [ 5 ] in Section V.C. 

B. Co-located EN and AP Deployment 

Next, we evaluate in Fig. |7] the performance of the proposed Algorithm [5] for co-located ENs and APs, 
where the value of Pr achieved by Algorithm [ 5 ] is plotted against the number of HAPs used (M). In 
particular, we compare its performance with that of LS (with M cluster centers as the initial searching 
points) and the CC placement method, i.e., the HAPs are placed at the M cluster centers. We can see that 
the proposed greedy algorithm in Algorithm [5] has the best performance among the methods considered. 
Nonetheless, the performance gaps over the LS and CC methods are relatively small compared to that 
in Fig. [51 An intuitive explanation is that the doubly-near-far phenomenon for co-located EN and AP 
renders the optimal HAPs placement to be around the cluster centers. By comparing Pig. [ 5 ] and Pig. [71 
we can see the evident performance advantage of using separated ENs and APs over HAPs. Por instance. 
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Fig. 8. Comparison of CPU time of the LS and Greedy methods for HAP placement optimization. 

the Pr achieved by 6 ENs and 6 APs in Fig. [5] is —O.lmW, while that achieved by 6 HAPs (equivalent 
to 6 ENs and 6 APs being co-located) is only around —O.lTmW. 

Although the greedy algorithm and the ES method perform closely, they differ significantly in the 
computational complexity. To better visualize the growth rate of complexity, we plot in Fig. [8] the 
normalized CPU time of the two methods, where each point on the figure is the normalized against 
the CPU time of the respective method when M = 4. Clearly, we can see that the complexity of the 
greedy algorithm increases almost linearly with M, where the CPU time increases approximately 6 times 
when M increases from 4 to 24. This matches our complexity analysis for Algorithm [5] in Section V.C 
that the complexity increases almost linearly in M when K is much larger. The ES method, however, 
has a much faster increase in complexity with M, where the CPU time increases by around 39 times 
when M increases from 4 to 24. Therefore, even in a large-size WPCN with large M, the computation 
time of the proposed greedy algorithm is still moderate, while this may be extremely high for the ES 
method, e.g., couple of minutes versus several hours for M = 50. 

C. Case Study: Separated Versus Co-located EN and AP 

Finally, we present a case study to compare the cost of node placement achieved by using either 
separated or co-located ENs and APs. Here, we consider a WPCN with 60 WDs uniformly placed within 
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Fig. 9. Comparison of the minimum deployment costs achieved by separated/colocated EN and AP placement optimization methods, where 
in (a) C 3 = 1.4 and in (b) cs = 1.05. 


the 24m x 24m box region, where the detailed placement is omitted due to the page limit. For the case 
of separated ENs and APs, we use Algorithm |4] to enumerate (M, N) pairs that can satisfy a given net 
energy harvesting performance requirement 7, and select the one with the minimum cost as the solution. 
For the case of HAPs, we use Algorithm [5] to find the minimum M that can satisfy the performance 
requirement 7. A point to notice is that, the obtained deployment solutions are sub-optimal to the min- 
cost deployment problems with either separated or co-located ENs and APs, i.e., the cost of the optimal 
deployment solution can be lower, because Algorithms IH and [5] are sub-optimal to solve problems (fTSl) 
and (fT5l) . respectively. More effort to further improve the solution performance is needed for future 
investigations. 

In Fig. Uta), we show the minimum deployment costs achieved by the two methods under different 
performance requirement 7. The number of nodes used by both methods are also marked in the figure. 
With {ci, Cg, C 3 } = {0.7,1,1.4}, we can see that using separated ENs and APs can achieve much lower 
deployment cost than co-located HAPs. In another occasion in Pig. |9}b), the two scheme achieve similar 
deployment cost when the cost of a HAP is decreased from C 3 = 1.4 to 1.05. We also observe that, by 
allowing the ENs and APs to be separately located, we need much less energy/information access points 
than they are co-located to achieve the same performance, thanks to the extra freedom in choosing both 
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the numbers and locations of energy/information access points. For instance, when 7 = 0, we need 24 
separated energy/information access points {M = 19 and N = 5), while 38 co-located energy/information 
access points by the 19 HAPs. However, we do not intend to claim that using separated ENs and APs is 
better than the co-located case. Rather, we show that a cost-effective deployment plan can be efficiently 
obtained using the proposed methods. In practice, the choice of using either separated or co-located 
ENs and APs depends on a joint consideration of the node deployment costs, the network size and the 
performance requirement. 


VIE Conclusions and Euture Work 

In this paper, we studied the node placement optimization problem in WPCN, which aims to minimize 
the node deployment cost while satisfying the energy harvesting and communication performance of 
WDs. Efficient algorithms were proposed to optimize the node placement of WPCNs with the ENs and 
information APs being either separated and co-located. In particular, when ENs and APs are separately 
located, simulation results showed that significant deployment cost can be saved by jointly optimizing 
the locations of ENs and APs, especially for large-size WPCNs with a large number of ENs and APs 
to be deployed. In the case of co-located ENs and APs (HAPs), however, the performance advantage of 
node placement optimization is not that significant, where we observed relatively small performance gap 
between optimized node placement solutions and that achieved by some simple heuristics, i.e., placing 
the HAPs at the cluster centers formed by the WDs. In practice, separated ENs and APs are more suitable 
for deploying WPCNs than co-located HAPs, because of the flexibility in choosing both the numbers and 
locations of ENs and APs. Nonetheless, because the optimal solution to the node placement problem has 
not been obtained in this paper, we may expect further improvement upon our proposed methods in the 
future, especially for the case of HAP placement optimization. 

Einally, we conclude with some interesting future work directions for the node placement problem in 
WPCN. Eirst, the models considered in this paper can be extended to more general setups. Using the uplink 
information transmission as example, we assumed in this paper that each WD has fixed association with 
a single AP. In practice, dynamic frequency allocation can be applied to enhance the spectral efficiency, 
where a WD can transmit to different APs, even multiple APs simultaneously, in different transmission 
blocks. Besides, instead of assuming each AP can serve infinite number of WDs, we can allow each 
AP to serve finite number of WDs. In addition, we may also consider the presence of uplink co-channel 
interference due to frequency reuse in WPCN. The extensions require adding corresponding constraints 


December 10, 2015 


DRAFT 


30 


or changing the expression of energy eonsumption model in the problem formulation of this paper. 
Seeond, it is interesting to eonsider the hybrid node deployment problem that uses both eo-loeated and 
separated ENs/APs. Third, it is praetieally important to eonsider the node plaeement problem with loeation 
constraints, e.g., some areas that may forbid the ENs/APs to be plaeed. In addition, the density of ENs 
may be eonstrained to satisfy eertain safety eonsideration on RE power radiation. 


Appendix 

Proof of the Convergence of Algorithm [3] 

Eet = !)••• ) -^1 denote the optimal solutions of (l20l) ealeulated from the /-th (/ > 1) set 

of assumptions made on the WD-AP assoeiations, denoted by k = 1- ■■ ,K. Eet denote the set 
of WDs to whieh the optimal solution vj^^’s eontradiet with the WD-AP assumptions (we eonsider only 
7 ^ 0, beeause otherwise the algorithm has reaehed its optimum), i.e.. 
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Aeeording to the proposed trial-and-error method, jk 
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Eet denote the minimum net energy harvesting rate among all the WDs given the updated WD-AP 


assoeiation and the current AP loeations v^^^’s, i.e.. 


^ _ Q _ a2,fe||v 0+1) - ) 

k=l-,K \ ’ 3k " / 


(35) 


We can see that beeause the update of WD-AP associations in (l34l) does not inerease the energy 

eonsumption rate of any WD achived by assuming jj^\ k = 1 ■ ■ ■ , K. Besides, j = 1, • • • , A^| 

is a feasible solution of (l20l) under the assoeiation assumption Therefore, the optimal solution 

|^(^+i) j = 1, • • ■ ,A^| ealeulated from the assoeiation assumption will lead to > 

f(*+i) > ^(0 In other words, the optimal objeetive of (l20l) is non-deereasing in eaeh trial-and-error update 
of WD-AP assoeiations. This, together with the faet that the optimal value of (l20l) is bounded, leads to 
the conclusion that the proposed trial-and-error method is convergent. 
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